Influence of high magnetic fields on superconducting transition of one-dimensional Nb 

and MoGe nanowires 
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The effects of strong magnetic field on superconducting Nb and MoGe nanowires with diameter 
~ fO nm have been studied. We have found that the Langer-Ambegaokar-McCumber-Halperin 
(LAMH) theory of thermally activated phase slips is applicable in a wide range of magnetic fields and 
describes well the temperature dependence of the wire resistance, over eleven orders of magnitude. 
The field dependence of the critical temperature, T c , extracted from the LAMH fits is in good 
quantitative agreement with the theory of pair-breaking perturbations that takes into account both 
spin and orbital contributions. The extracted spin-orbit scattering time agrees with an estimate 
t so ~ T^hc/Ze 2 ) 4 , where r is the elastic scattering time and Z is the atomic number. 

PACS numbers: 74.78.Na, 74.25. Fy, 74.25.Ha, 74.40.+k 



The problem of superconductivity in one- dimensional 
(ID) systems attracts much attentions since it involves 
such fundamental phenomena as macroscopic quantum 
tunnelling, quantumphase transitions and environmen- 
tal effects]! E H II IE It is expected that a strong 
magnetic field can be used to control these phenomena. 
Indeed, the microscopic theory predicts that a magnetic 
field, acting on a superconducting condensate, lifts the 
time reversal symmetry of the spin and orbital states of 
paired electrons and suppresses the critical temperature, 
T c 0, H . A strong enough field destroys superconductiv- 
ity. The magnetic field pair-breaking effects were studied 
in de pth in two and zero- dimensional systems, i.e. thin 
films | id] and nanograins [TJ]. However, an experimental 
verification of the pair breaking effects in ID supercon- 
ductors is long overdue. 

A distinct feature of ID superconductors is the absence 
of the phase coherence. Due to fluctuations the ampli- 
tude of the order parameter has a finite probability to 
reach zero at some point along the wire, allowing the 
phase of the order parameter to slip by 2-7T |l2( . The the- 
ory of thermally activated phase slips (TAPS) was devel- 
oped by Langer, Ambegaokar, McCumber and Halperin 
(LAMH). However the effect of the magnetic field on the 
phase slippage process is not established. It is also un- 
known whether the magnetic field can change the relative 
contributions of quantum and thermally activated phase 
slips in thin wires 0, 0] • 

In this Letter we study the effects of the magnetic field 
on the phase slippage rate and the critical temperature 
of thin wires. It is found that the LAMH provides a good 
description for ID superconductors in magnetic fields up 
to 11 T. The dependence of the critical temperature on 
the magnetic field, T C (B) agrees well with the theory of 
pair-breaking perturbations that takes into account both 
spin and orbital contributions [8|, [9| • This is our main 
result. No significant contribution of quantum phase slips 
has been detected in the studied samples. 

The samples were made by sputter- coating of sus- 



pended fluorinated carbon nanotubes with Nb or 
Mo7gGe2i. Transport measurements wereperformed in a 
He-3 cryostat, as described previously jj,y,|5{. The mag- 
netic field was oriented perpendicular to the wire and 
parallel to the thin film electrodes connected in series 
with the wire. 

A series of resistance versus temperature R(T) curves 
measured at different magnetic fields is shown in Fig. 1. 
For each curve, a resistance drop at higher temperature 
corresponds to a superconducting transition in the film 
electrodes. The resistance value immediately below the 
drop is taken as the wire normal resistance Rn- The 
second resistive transition corresponds to development of 
superconductivity in the wire. With increasing magnetic 
field, both transitions shift to lower temperatures. 

To analyze these data we employ the LAMH expression 




FIG. 1: Resistance versus temperature dependence for a 
Nb nanowire (the thickness is 8 nm, the normal resistance 
is R N = 700 Q,, the length is L=120 nm). Each R(T) curve is 
measured in a fixed magnetic field; some fields are indicated. 
Solid lines show the fits to the LAMH theory. 
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where AF = (8 v / 2/3)(if c 2 /87r)^ is the energy bar- 
rier, CI = (L I '^(AF/kBT) 1 / 2 '(I /t G l) is the attempt fre- 
quency, and tql = [7f/i/8fcs(r c — T)] is the Ginzburg- 
Landau (GL) relaxation time, L is the length of the wire, 
A its cross-sectional area, and £ is the GL coherence 
length. Following Ref. 3] we express the energy bar- 
rier as 



AF(T) « 0.8Z[L/£(0)]{R q /R N )k B T c {l - T/T c 
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where R q = h/4e 2 = 6.45 kQ is the resistance quan- 
tum for Cooper pairs and £(0) is defined by = 
£(0)(1 — T/T c ) -1 / 2 . Taking into account the contribu- 
tion of quasiparticles, the total resistance is given as 
R = (R^ 1 + R~lamh)~ 1 '■ I n our fitting procedure we 
use two adjustable parameters, T c and £(0). 

The LAMH fits are shown as solid lines in Fig. 1, for 
various magnetic fields. Although the LAMH theory is 
derived for B=0, we find that the resistance agrees with 
the LAMH fits very well, for both Nb and MoGe samples, 
even in high fields up to ~11 T. 

Our extension of the LAMH theory to high magnetic 
fields requires an explanation: For a phase slip to occur 
in a wire, the system needs to overcome an energy bar- 
rier that is a product of the condensation energy density 
(H 2 /8n) in a volume of a phase slip A£. It can be shown 
within the GL theory that in magnetic field the conden- 
sation energy density goes to zero as (1 — T/T C (B)) 2 , 
where T C (B) is the field-dependent critical temperature. 
The coherence length varies as £(0, B)(l — T/T C (S)) -1 / 2 
and diverges at T C (B). Because the temperature depen- 
dence of both the condensation energy and the coherence 
length has the same form as in zero field, we expect Eq. 
1 and Eq. 2 to be applicable in magnetic fields also. The 
observed agreement with the experiment suggests that 
the mechanism of the phase slippage in ID wires is not 
changed by magnetic field. 

To test the LAMH theory in a nonlinear regime, we 
performed measurements of the volt age- current V(I) de- 
pendencies at high bias currents for some MoGe wires 
(Fig. 2a). We observe that at high currents the wires 
undergo a transition into resistive state. This transition 
is smooth at high temperatures and it is jumpwise, with 
some hysteresis, at low temperatures 0, 0] (Fig. 2a). 
A closer inspection of the data revealed that even at low 
temperatures there exists a small non-linear voltage vari- 
ation at currents slightly lower than the switching cur- 
rent (Fig. 2b). The V(I) curves remain qualitatively 
unchanged in magnetic field. 

In order to investigate in more detail the resistive tails 
observed slightly below the switching current we mea- 
sured the differential resistance dV/dl versus bias cur- 
rent / as shown in Fig. 3. It is clear that at all tem- 
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FIG. 2: (a)Voltage versus current dependence at indicated 
temperatures for a MoGe nanowire (thickness 7 nm, Rn = 4.3 
kfi, L=190 nm) at B=0. The curves are vertically shifted by 
1 mV for clarity, (b) The same data at a magnified scale. 



peratures the dV/ dl versus / data follow exponential de- 
pendence, which is expected from the LAMH expression 
dV/dl = R(T) cosh (J/ Jo), where R(T) is the zero-bias 
resistance given by Eq. 1 and Iq = AekT/h. We fit 
the data with the above expression and extract two ad- 
justable parameters R(T) and Iq. Experimental values of 
Iq are close to the theoretical value Iq = AekT/h (Fig. 4, 
insert). We speculate that the observed small upward de- 
viation of Iq might be due to the fact that some fraction 
of the bias current is carried by nonequilibrium quasi- 
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FIG. 3: Differential resistance versus bias current at indicated 
temperatures for the MoGe nanowire (thickness 7 nm, Rn = 
3.9 kfi, L=I50 nm) in a fixed magnetic field B=9 T. Solid lines 
are fits to the LAMH expression dV/dl = R(T) cosh(/// ) 
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FIG. 4: Resistance versus temperature for the MoGe nanowire 
(same as in Fig. 3) in magnetic fields and 9 T. Open circles 
represent zero-bias-current measurements and black squares 
indicate the resistance values obtained from the fit of the non- 
linear portion of dV/dl curves (Fig. 3). Solid lines are the 
fits to the LAMH theory (Eq.l). Extracted fitting parameters 
T c and £(0) are indicated. The insert shows the experimental 
dependence of the parameter Iq on temperature (solid and 
open symbols) and the theoretical value Iq = 4ekT/h (solid 
line). 



particles. Such quasiparticles are generated by the phase 
slips and are not taken into account within the LAMH. 
A detailed theoretical analysis is needed for further un- 
derstanding of the Iq(T) behavior. 

In Fig. 4 we superimpose the zero-bias resistance data 
and the resistance data obtained from the fits to the non- 
linear portion of dV/dl vs. I curves. These two sets of 
data appear mutually consistent. The LAMH fit (solid 
lines in Fig. 4) gives an excellent description of all re- 
sistance data in a range of eleven orders of magnitude. 
We also find a good agreement with the LAMH for data 
taken in magnetic field B=9 T (Fig. 4). It is therefore 
concluded that the resistance in studied nanowires is de- 
termined by thermally activated phase slips even in high 
magnetic fields. 

In Fig. 5 we plot the extracted parameters T C (B) and 
£(0, B) versus magnetic field. The coherence length in- 
creases very slowly in low fields and starts to grow more 
rapidly at B > 8 T. We consider now the critical tem- 
perature T C (B), which is the main focus of this work. 
Although the initial decrease of T c agrees with the varia- 
tion predicted by the GL theory (Fig. 5a), it is clear that 
the GL phenomenology is not sufficient to account for the 
observed T C (B) dependence. One possible reason is that 
the GL theory does not take into account the spin pair- 
breaking effect and the coefficients of the GL theory can 
change in the high field regime. Thus we have to use the 
exact theory of pair-breaking 8] . In both Nb and MoGe 



FIG. 5: Adjustable parameters used in Eq.f for the LAMH 
fits of the type shown in Fig.l. (a) The critical temperature 
versus magnetic field T C (B) is shown for the Nb wire (same as 
in Fig.l) and for a MoGe wire, (b) The GL coherence length 
versus magnetic field. The parameters of the MoGe wire are: 
thickness is 9 nm, Rn = 7.5 kf2, L=460 nm. The theoretical 
Ginzburg-Landau dependence is indicated by dashed lines. 



nanowires superconductivity persists to magnetic fields 
that are larger than the paramagnetic limit, _B P [T]=1.84 
T c [K] 9], (10.6 T for Nb and 8.1 T for MoGe). Since the 
superconductivity is not fully suppressed at such fields, 
we conclude that the effect of magnetic field on the spin 
part of the Cooper pair is reduced by the spin-orbit scat- 
tering, as expected for materials with high atomic num- 
bers. If the spin-orbit scattering is sufficiently strong, 
the transition into the normal sate is continuous. This 
allows us to take into account both the spin and the or- 
bital effects on T c in an implicit relation of the theory of 
pair-breaking perturbations 

, T C (B) fl\ , (I a + a, \ 

where ijj(z) is the digamma function and a Q = 
2De 2 (A 2 )/hc 2 and a s « r so e 2 hB 2 /m 2 c 2 are the or- 
bital and the spin pair-breaker strength parameters, D 
is the diffusion coefficient and r so is the spin-orbit scat- 
tering time. To find the vector potential averaged over 
the cross-sectional of a wire, (A 2 ), we use the ex- 
pression for a cylinder in perpendicular magnetic field, 
(A 2 ) = B 2 d 2 /16, where d is the diameter of the wire. 

By solving Eq. 3 numerically the theoretical depen- 
dence of the normalized critical temperature versus the 
normalized magnetic field is obtained. In Fig. 6 this 
dependence is compared to the experimental values of 
T C (B) /T c (0) (same data as in Fig. 5a). Since both pair- 
breakers, a a and a 8 , have a quadratic dependence on £?, 
only one adjustable parameter is required, which is the 
critical field of the wire at zero temperature, B wc . The 
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H3, £=1.2 nm 0), and so t so ~ 2.4 x 10 -13 s, again 
in agreement with the experimental value. Thus we con- 
clude that the pair-breaking theory combining spin and 
orbital contributions gives an accurate quantitative pre- 
diction for the suppression of the critical temperature of 
homogeneous ID superconducting wires. 

This work is supported by NSF CAREER Grant No. 
DMR 01-34770, by A.RSloan Foundation and by a DOE 
grants DEFG02-91ER45439 and DEFG02-96ER45439. 



FIG. 6: Normalized critical temperature versus normalized 
magnetic field for the Nb and MoGe wires. Critical magnetic 
fields (B cw ) for both wires are indicated. The solid line is a 
fit to the theory of pair- breaking perturbation (Eq. 4). 



best fits are obtained by choosing B wc =11.7 T for Nb 
and 16.5 T for MoGe (Fig. 6). 

To compare orbital and spin pair-breaking contribu- 
tions we introduce an orbital critical field (B co ) and 
a spin critical field (B cs ), defined as fields needed to 
suppress superconductivity if only one of the two pair- 
breaking mechanisms is present. These fields are ob- 
tained from equations 2a a = 2a s = 1.76/cT c (0) 9j. The 
total critical field is then B wc = (B~ 2 + S" 2 )" 1 / 2 . Us- 
ing BCS and GL relations the orbital critical field can 
be written as B co = 0.53<&rj/^£(0), where $o is the flux 
quantum. The value of the coherence length £(0) is 
known from the LAMH fit in zero field (7.1 nm for Nb 
and 7.7 nm for MoGe). Taking as d the thickness of 
the deposited material (corrected for oxidation for the 
MoGe wire 5J) we estimate B co m 22 T for both wires. 
From B co and experimental total critical field B wc we 
determine the spin critical fields, B cs « 14 T for Nb and 
B cs « 25 T for MoGe. Thus the spin and orbital pair- 
breakers have comparable strengths. 

From the B cs we estimate the spin-orbit scattering 
time as r so = 2.3 ± 0.5 x 10 -13 s for the Nb wire and 
t so = 5±3 x 10" 14 s for the MoGe wire. The MoGe result 
is considerably different from the value r so ~ 1.3 x 10 -12 
s obtained for thin MoGe films from weak localization 
measurements [u|. With such r so value the supercon- 
ductivity in MoGe wire would be completely suppressed 
at B cs ~ 5 T, contrary to our observation. On the other 
hand, we can use the formula t so ~ r(Hc/ 'Ze 2 ) A given in 
Ref. [u| . With elastic scattering time r ~ 6 x 10~ 16 s [u| 
this gives a shorter spin-orbit scattering time 8 x 10 -14 
s, that agrees with our result. The latter estimate also 
works well for the Nb nanowire. Here we have elastic 
time r = £/v F ~ 1.9 x 10" 15 s (v F = 0.62 x 10" 8 cm/sec 
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